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We apply the escape-rate formalism to compute the shear viscosity in terms of the chaotic proper- 
ties of the underlying microscopic dynamics. A first passage problem is set up for the escape of the 
Helfand moment associated with viscosity out of an interval delimited by absorbing boundaries. At 
the microscopic level of description, the absorbing boundaries generate a fractal repeller. The fractal 
dimensions of this repeller are directly related to the shear viscosity and the Lyapunov exponent, 
which allows us to compute its values. We apply this method to the Bunimovich-Spohn minimal 
model of viscosity which is composed of two hard disks in elastic collision on a torus. These values 
are in excellent agreement with the values obtained by other methods such as the Green-Kubo and 
Einstein-Helfand formulas. 

PACS numbers: 05.45.Ac; 05.45.Df; 05.45.Jn; 05.60.-k] 

I. INTRODUCTION 

The link between the irreversible phenomena governed by macroscopic equations such as the Navier-Stokes equations 
and the microscopic reversible dynamics of the atoms and molecules is a fundamental problem. In this context, it 
has been shown that typical many-body systems of interacting particles present a chaotic dynamics P, 0, 0) • This 
microscopic chaos develops a sensitivity to initial conditions over a time scale of the order of the intercollisional time 
of the atoms and molecules. This sensitivity to initial conditions is characterized by positive Lyapunov exponents 
which are the rates of exponential separation between some reference and perturbed trajectories of the system. The 
sensitivity to initial conditions results into a huge dynamical randomness characterized by a positive Kolmogorov-Sinai 
(KS) entropy per unit time given by the sum of positive Lyapunov exponents if the system is at equilibrium: 

\i>0 

an identity known as Pesin's equality 0, 0] . This dynamical chaos provides an efficient mechanism of randomization 
of the different observable quantities such as the microscopic currents associated with the transport properties. 

For a Hamiltonian-type microscopic dynamics, a connection can be established between dynamical chaos and the 
transport properties thanks to the escape-rate formalism 0, 0j ^i ^| • In this formalism, the gap between the 
kinetic time scale of chaotic properties and the hydrodynamic time scale of transport properties is bridged by linking 
the transport coefficients to differences between chaotic properties. In this formalism, the following formula has been 
derived for viscosity ^] 

-/iKsj , (2) 

where the difference between the sum of positive Lyapunov exponents and the KS entropy is nonvanishing because 
the system is here under nonequilibrium conditions. These nonequilibrium conditions select the trajectories of the 
many-body system which do not escape out of a phase-space region specific to the transport property of interest. This 
phase-space region is defined by requiring that the Helfand moment associated with the transport property remains 
bounded in an interval of extension x- In the limit x ~^ oo, the nonequilibrium condition is progressively relaxed and 
the sum of positive Lyapunov exponents as well as the KS entropy approach their equilibrium value satisfying Pesin's 
equality Under nonequilibrium conditions, Pesin's equality is not satisfied and the difference gives the rate of 
escape of trajectories out of the aforementioned phase-space region 0^3- This region contains a fractal repeller J^^ 
composed of trajectories which escape neither in future nor in past. The escape rate is characteristic of this fractal 
repeller and is related to the transport coefficient, leading to the formula ||2J). The escape-rate formalism has already 
been applied to the transport property of diffusion as well as to reaction-diffusion processes [T3 . ITsI . 
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The purpose of the present paper is to apply the escape-rate formahsm to viscosity. The syst em we use as a vehicle 
of our study is a minimal model of viscosity previously analyzed by Bunimovich and Spohn [la | . The minimal models 
of transport are of special interest because they are the simplest possible models already presenting a positive and 
finite transport coefficient. It is known that a minimal model should contain only one particle for diffusion, two 
particles for viscosity, and three particles for heat conduction (16.j . For viscosity, we therefore consider here the model 
composed of two hard disks moving on a torus and undergoing elastic collisions. In a previous paper |l7j . we have 
described this model and some of its properties for a hexagonal and a square geometries. Our aim is here to compute 
viscosity thanks to the escape-rate formalism and to show the equivalence with the results of the Green-Kubo and 
Einstein-Helfand formulas already obtained in the previous paper . 

In the present paper, we compute the viscosity by using the chaotic and fractal properties of the repeller. We use 
a variant of Eq. |(2Jl in which the difference between positive Lyapunov exponents and KS entropy is given in terms 
of the Lyapunov exponents and the partial fractal dimensions of the fractal repeller. Indeed, the fractality of the 
repeller is a corollary of its chaoticity so that its fractal dimensions are related to its KS entropy. Accordingly, the 
knowledge of the partial dimensions allow us to evaluate the KS entropy. 

The paper is organized as follows. In Sec. we develop the escape-rate formalism for shear viscosity. In Sec. IIIII 
we present the fractal repeller of viscosity in the two-disk model, which we compare with the fractal repeller of diffusion 
in the Lorentz gas. In this way, we show that a specific fractal repeller is associated with each transport property. 
Finally, the chaotic and fractal properties are described in Sec. IIVI where we compute the viscosity coefficient from 
the positive Lyapunov exponent and the fractal dimension of the repeller. Conclusions are drawn in Sec. 



II. VISCOSITY IN THE ESCAPE-RATE FORMALISM 



A. Helfand moment for viscosity 



In the previous paper p^ . we have shown that the shear viscosity coefficient can be computed with the Einstein- 
Helfand formula 



where the Helfand moment is defined by 



Vyx,yx = lim — {[Gyx{t) - Gy:r(0)]^ 
t^oo Zt 



N 
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a=l 



(3) 



(4) 



Gyx = \J Y — \l T7 

for a system of N particles of position = (xq, ...) and momentum — (pn ^iPav, ■■■) moving in a domain 
delimited by periodic boundary conditions. As explained in the previous paper [131 , the particles must satisfy the 
minimum image convention, which requires the presence of the extra terms in the Helfand moment (0J) involving the 
jumps Ax'^^ of the particles to fulfil the periodic boundary conditions, are the times of the jumps, pii^y = Payitg) 
is the momentum at the time tg of the jump. In Eq. (Q, the average (•) is taken in the equilibrium microcanonical 
ensemble for which 



N 

kBT{N-l) 



(5) 



B. First passage problem for viscosity 



The central object of the escape-rate formalism is the fractal repeller which is composed of the phase-space trajec- 
tories for which the Helfand moment fluctuates forever (in the future and the past) within some interval: 

- I < G,. < -t-| . (6) 

These trajectories are exceptional because the Helfand moment escape out of this interval for almost all the trajectories. 
Therefore, the repeller has a vanishing probability measure in the phase space albeit it is typically composed of a 
nonenumerable set of trajectories. The repeller thus typically forms a fractal in the phase space |9l[l0|. 

We set up a first passage problem of the Helfand moment by introducing absorbing boundaries at Gyx = These 
absorbing boundaries in the space of variations of the Helfand moment correspond to equivalent absorbing boundaries 
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in the phase space of the system. In the phase space, the absorbing boundaries dehmit a domain which contains the 
fractal repeller. We consider a statistical ensemble of initial conditions taken inside this domain and we run their 
trajectories. When a trajectory reaches the absorbing boundaries it escape out of the domain and it thus removed 
out of the statistical ensemble. 

Under the forward time evolution, the remaining trajectories belong to the stable manifolds of the repeller. Under 
the backward time evolution, the remaining trajectories belong to the unstable manifolds of the repeller. Under 
both the forward and backward time evolution, the remaining trajectories belong to the repeller itself which is the 
intersection of its stable and unstable manifolds 9]. For a typical chaotic dynamics, almost all trajectories escape out 
of the domain after some time so that the repeller as well as its stable or unstable manifolds are fractal objects. 

These fractals can be generated by allowing the escape of trajectories over a long but finite time interval. Over a 
finite time, there remain a sizable set of trajectories, which progressively reduces to the fractal as the time interval 
becomes longer and longer. The number of trajectories in the set (or statistical ensemble) decays with time. Typically, 
the decay is exponential and characterized by the so-called escape rate 7. 

The escape rate 7 can be evaluated by solving the first passage problem of the Helfand moment by introducing 
absorbing boundaries at Gyx = ±f • Indeed, the Einstein-Helfand equation (jHJ shows that the Helfand moment 
performs a diffusive-like random walk. Accordingly, the Helfand moment can be considered as a random variable 
g — Gyx for which the probability density p((;) obeys to a diffusion-type equation Q: 

where the role of the diffusion coefficient is played by the shear viscosity Q itself. At the absorbing boundaries, the 
probability density must satisfy the absorbing boundary conditions: 

The solution of the diffusion-type equation iQ with the boundary conditions © is given by 



i) = ^ c„ exp(-7„i) sin 

with 



X V 2 



(9) 
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and where the coefficient c„ depend on the initial probability density. The number Af{t) of trajectories remaining 
between the absorbing boundaries at the current time t is related to the probability density by: 

Mit) = M) / Pig, t) dg ~ Afo exp(-7it) . (11) 



After a long time, the escape is dominated by the smallest decay rate, 71, which can therefore be identified with the 
escape rate 7. In this way, we obtain the escape rate as a function of x- 
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7 = 7i=77(^-j . (12) 

This result is obtained by using the diffusion-type equation Q which is expected to hold over spatial distances larger 
than the mean free path of the particles. Therefore, the parameter x must be sufficiently large so that the Helfand 
moment is in a diffusion regime and Eq. Q holds. 

The shear viscosity coefficient can thus be obtained from the escape rate which depends on the parameter x of 
separation between the absorbing boundaries as 



V = lim (-) 7(x) ■ (13) 



In the following we call Eq. (|13|l the escape-transport formula. 
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C. The chaos-transport formula 



At the microscopic level of description, the escape rate is controlled by the fractal repeller which is composed 
of all the trajectories which satisfy the condition (|SJ| under forward and backward time evolution. The repeller is 
the support of a natural invariant probability measure. This invariant measure is natural because it is generated 
by the dynamics and can be approximated by a statistics based on the trajectories remaining within the absorbing 
boundaries after a long but finite time. The dynamics on the fractal repeller is characterized by positive Lyapunov 
exponents and a KS entropy, both evaluated with respect to the natural invariant measure of the repeller. If the 
dynamics is unstable some Lyapunov exponents are positive. If the dynamics is chaotic the KS entropy is positive. 
On a repeller, the sum of positive Lyapunov exponents differs from the KS entropy and the difference gives the escape 
rate 7(x) of the repeller T^: 



7(X) = ( E^'-^Ks) 

\A,>0 / 



(14) 



If we combine this result from dynamical systems theory with Eq. we obtain the chaos-transport formula ^ 

for viscosity as originally derived by Dorfman and Gaspard j^. 

An equivalent formula can be obtained which involves the partial fractal dimensions of the repeller instead of the 
KS entropy. Indeed, the fractal character of the repeller is a direct consequence of the escape of trajectories so that 
the KS entropy is no longer equal to the sum of Lyapunov exponents but to 

/iKS = E Ai , (15) 

Ai>0 

where the coefficients are the partial information dimensions of the repeller associated with each unstable direction 
of corresponding Lyapunov exponent '3| . These partial dimensions satisfy 

< < 1 , (16) 

so that the KS entropy is in general smaller than the sum of positive Lyapunov exponents. Accordingly, the escape 
rate can be expressed as 

\Ai>0 / 

in terms of the partial codimcnsions defined as 

c. = 1 - d^ . (18) 
Combining with Eq. I|13|l . the shear viscosity is given by 



\A,>0 / 



(19) 



In the limit x ^ the Lyapunov exponents reach their equilibrium values Ai^oq, while the codimcnsions vanish 
typically as q ^ if transport is normal. If we introduce the quantities 

a, = lim f-)%^ , (20) 

Eq. (|19|) provides a decomposition of the viscosity coefficient on the spectrum of Lyapunov exponents such as 

Ai,oq>0 

Typically, the escape is most important in the most unstable direction corresponding to the maximum Lyapunov 
exponent Ai. Therefore, the repeller is more depleted in the most unstable direction and the corresponding partial 
dimension d\ is lower than the further ones. This reasoning suggests that a typical behavior is 

(ci Ai)^^ , (22) 
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for X ^ oo if the maximum Lyapunov exponent Ai is well defined. 

This is precisely the case in two-degree-of-freedom systems such as the two-disk model where the chaos-transport 
formula reduces to 

lim (^)'(ciA)^ , (23) 

where A is the unique positive mean Lyapunov exponent and ci the corresponding codimension which should be 
understood as the partial information codimension of the unstable manifolds of the fractal repeller given in terms of 
the partial information dimension by Young's formula [T^] 

CI = 1 - di = 1 - . (24) 

It is known that the partial information dimension of the repeller is well approximated by the partial Hausdorff 
dimension if the escape rate is small enough and if Ruelle's topological pressure does not present a discontinuity. 
This last condition is fulfiled if the system does not undergo a dynamical phase transition, which is the case in the 
finite-horizon regimes of Sinai's billiard which controls the dynamics of both the Lorentz gas and the two-disk model 
|13|. Under these conditions, we can replace the partial information codimension by the partial Hausdorff codimension 
in the chaos-transport formula and obtain the viscosity as 

V= lim (^)'(cH A)^ . (25) 

In the limit x ^ co? the Lyapunov exponent converges to its equilibrium value so that we can also write: 

77 = Acq lim (-) ch(x) , (26) 

under similar conditions as Eq. (|25|l . 



III. FRACTAL REPELLER 



In this section, our purpose is to display the fractal repeller associated with viscosity in the two-disk model and 
to compare it with the fractal repeller of diffusion in the Lorentz gas in order to show that they are different and 
therefore specific to each transport property. 



A. Shear viscosity in the two-disk model 

The two-disk model has been studied by Bunimovich and Spohn who proved thanks to a central limit theorem that 
viscosity is positive and finite in this minimal model [l^ . In the previous paper JJJ , we have considered the two-disk 
model in the hexagonal geometry which we shall use in the following. We showed in Ref. [l3 that the dynamics 
reduces to a Sinai billiard in the center-of-mass frame and that the Helfand moment l@J with A'^ = 2 is then given by 



Gyx{t) 



ts) 



(27) 



where {x, y) are the coordinates of the relative position of both disks and {px,Py) the canonically conjugated relative 
momentum. The jumps happen when the trajectory of Sinai's billiard crosses the hexagonal boundary. If the trajectory 
crosses the side of label uj the trajectory is reinjected at the opposite side so that the jump in position is given by the 
lattice vector Ar*^*^ — — cif-* corresponding to the side uj. 

A fractal repeller is defined by considering all the trajectories such that their Helfand moment satisfies the conditions 



(28) 



where the parameter x should be large enough. The stable manifolds of the fractal repeller can be visualized by 
plotting the initial conditions of trajectories satisfying the conditions H28|) over a long time interval extending forward 
in time. These initial conditions are taken on the disk of Sinai's billiard. The initial conditions are specified by the 
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FIG. 1: Initial conditions of the particules in the Sinai billiard. 



angle 9 of the initial position and the angle <j} that the initial velocity makes with a vector which is normal to the disk 
at the initial position (see Fig. The initial conditions are plotted in the Birkhoff coordinates (0,sin0). 

Figm'e|21depicts such the fractal composed of the stable manifolds of the repeller for viscosity in the two-disk model. 
We provide evidence that the set is fractal by zooming successively on it in Figs. O and 01 where the self-similarity of 
the repeller clearly appears. 

Let us take a section across the repeller in Fig. |21 at Oq — it /A. Taking the escape time of the corresponding 
trajectory, we have obtained the escape-time function depicted in Fig. [S] The time for the trajectory to escape 
out of the phase-space region corresponding to the interval (|28|) is infinite if the trajectory belongs to the stable 
manifold of a trajectory of the repeller. Indeed, this trajectory is then asymptotic to a trajectory which does not 
escape. Accordingly, the escape-time function has vertical asymptotes on the stable manifolds of the repeller. Since 
the repeller is fractal the vertical asymptotes are not enumerable, which explains the behavior in Fig. |S| 







FIG. 2: Fractal repeller associated with viscosity in the hexagonal geometry with absorbing boundaries at x = 2.70. The 
density is n = {2/V) = 0.45. 
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FIG. 3: Enlarging of the domain in to the square in Fig. |5| 



B. Diffusion in the Lorentz gas 

Diffusion of a tracer particle in the hard-disk periodic Lorentz gas has been studied with the escape-rate formalism 
in Ref. In this Lorentz gas, the tracer particle undergoes elastic collisions on hard disks forming a triangular 

lattice. In a unit cell of the lattice, the dynamics also reduces to Sinai's billiard. The energy of the tracer particle 
is conserved as well as the phase-space volumes. Sinai and Bunimovich have proved that the dynamics is ergodic 
and mixing and that the diffusion coefficient is positive and finite in the finite- horizon regime |20|. For diffusion, 
the associated Helfand moment is simply given by one of the coordinates {x,y) of position of the tracer particle j^]- 
An escape process is associated with diffusion by setting up a problem of first passage of the tracer particle at some 
absorbing boundaries. If we consider the x-coordinates, the tracer particle does not escape as long as the following 
condition is satisfied: 

-^<x<+-. (29) 

The absorbing boundary conditions are therefore defined at a; = ±-f . With these absorbing boundaries, the system 
is called an open Lorentz gas [l^ . 

The trajectories trapped within the interval H29(l form a fractal repeller as shown in Ref. In order to compare 
with the fractal repeller of viscosity, we can plot the fractal repeller of diffusion in a similar way as here above for 
viscosity. 

Here again, we plot all the initial conditions of trajectories remaining within the interval (|29|l over a long forward 
time interval. These initial conditions are plotted in the same Birkhoff coordinates (6', sin cj)) of a disk around the 
coordinate a; ~ in the Lorentz gas. The set of the selected initial conditions approximate the stable manifolds of the 
fractal repeller. We zoom successively on this fractal in Figs. [7|and|Sl which provides evidence of its self-similarity. 
As a consequence, the repeller is also fractal. The fractal dimension of the repeller is related to the diffusion coefficient 
of the Lorentz gas and its Lyapunov exponent, as shown in Refs. 0i0|- 
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FIG. 4: Enlarging of the domain into the square in Fig. |3 



C. Comparison between diffusion and viscosity 

We observe that the two fractal repellers associated respectively with diffusion (see Fig. and viscosity (see Fig. 
121) are different. Indeed, although the global structure is similar, the trajectories belonging to the different repellers 
are not the same. 

To convince us of this difference, we take some examples of trajectories. In Fig. we have a periodic trajectory 
bouncing between two disks in the billiard. This trajectory belongs to the repeller associated with diffusion since the 
position X is bounded and satisfied (|29|l . However, the viscosity Helfand moment of this trajectory does not satisfy the 
condition H28() so that it does not belong to the repeller of shear viscosity. With Eq. (|27|l . we see that, in one direction, 
both Aa;'-'*'' and p^'^ are positive. Therefore, the contribution at this passage is positive for the Helfand moment. In 

the other direction, both Aa;^**^ and p^^ are negative but the product Ax^^^py'^ is also positive. Accordingly, the 
Helfand moment increases forever on this trajectory which, therefore, does not belong to the repeller associated with 
shear viscosity. 

On the other hand, we can observe the opposite case. Figure [TUI depicts an example of trajectory escaping from the 
interval (|29|l although its Helfand moment of viscosity remains between the absorbing boundary conditions (|28|l . 
The repellers associated with different transport properties are therefore different. 



D. Escape rate and viscosity 

In this subsection, we show that the shear viscosity can be obtained from the escape rate of the repeller by using 
the escape-transport formula . We consider a sequence of repellers with larger and larger values of the parameter 
X- The escape rate 7(x) is numerically evaluated for each repeller by computing the decay of the number M{t) of 
trajectories still within the interval H28|l at current time and by extracting the escape rate from the exponential decay. 
The escape rate is observed to behave as 7(x) ^ and the shear viscosity coefficient is then obtained with Eq. 

US; 

Figures 1111 and 1121 depict the viscosity directly computed from the escape rate and compared with the values 
obtained by the Einstein-Helfand formula in Ref. (13] i respectively in the hexagonal and square geometries. As in the 
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FIG. 6: Fractal repeller associated with diffusion in the hexagonal geometry with absorbing boundaries at _R = 4. The density 
of hard disks is n = 0.45. 



previous paper |17|. we consider reduced viscosities defined by 



7?*, = -JM=^ . (30) 



We observe in Figs. 1111 and [T^ the excellent agreement between both methods. 
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FIG. 7; Enlarging of the domain into the square in Fig. |S| 

IV. VISCOSITY FROM THE CHAOTIC AND FRACTAL PROPERTIES OF THE REPELLER 

In this section, we compute the shear viscosity coefficient in terms of the chaotic and fractal properties of the repeller 
by using the chaos-transport formula (|25|) which related the viscosity to the Lyapunov exponent and the Hausdorff 
codimension of the repeller of viscosity. 

A. Lyapunov exponent 

In Sinai's billiard which controls the reduced dynamics of the two-disk model, the elastic collisions between the 
disks are defocusing. This induces a dynamical instability of the trajectories which is characterized by the Lyapunov 
exponents. These exponents are the rates of exponential separations between a reference orbit and infinitesimally close 
orbits. Since the dynamics of Sinai's billiard is symplectic and volume-preserving in the four-dimensional phase space, 
the Lyapunov exponents spectrum is (-f A, 0, 0, —A) so that their sum is vanishing. One of the Lyapunov exponents 
vanishes because of the absence of exponential separation in the direction of the flow. Another one corresponding to 
the direction perpendicular to the energy shell equals zero because of energy conservation. 

There exists a method to calculate the positive Lyapunov exponent by considering the motion of a front of particles 
accompanying the reference particle and issued from the same initial position but with different initial velocities |2lj | . 
Because the dynamics is defocusing, this front is expanding. Locally on the reference orbit Ft the front (called the 
unstable horocycle) is characterized by a curvature Kui^t) or, equivalently, by its radius of curvature 1/K„(rt). Thanks 
to this method explained in detail in Ref . [Tol Il3j| , we have computed the positive Lyapunov exponent as a function of 
the density of the system (in the hexagonal and square geometries) . The equilibrium values of the Lyapunov exponent 
are obtained by running a trajectory in Sinai's billiard without absorbing boundaries and by averaging over a long 
time interval. The resulting numerical values are depicted in Figs. 1131 and [T^ 

In the chaos-transport formula (|25|) . the Lyapunov exponent has to be evaluated for the trajectories belonging to the 
fractal repeller. The statistical average is here carried out for the natural invariant probability measure concentrated 
on the fractal repeller. This invariant measure defines a nonequilibrium state for the motion. As aforementioned, the 
natural invariant measure is generated by the dynamics itself. Accordingly, the Lyapunov exponent is numerically 
computed by averaging over a statistical ensemble of trajectories which has not yet escaped after a long but finite 
time. This ensemble can be as large as wished by increasing the number of initial conditions. In this way, we can 
calculate the nonequilibrium values of the Lyapunov exponent. 
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FIG. 8; Enlarging of the domain into the square in Fig. Q 




FIG. 9: Periodic trajectory belonging to the fractal repeller associated with diffusion but not to the one associated with 
viscosity. 



In Table I, we present a comparison between the equilibrium Lyapunov exponent Acq without absorbing boundary 
conditions (as depicted in Fig. ^| and El and the nonequihbrium Lyapunov exponent Aneq(x) evaluated over a 
nonequilibrium measure which has the fractal repeller as support. The difference between these exponents is small 
and of the order of the escape rate, in agreement with the results of Ref. for the disordered Lorentz gas. 



Table I. Values of the characteristic quantities of chaos for different densities n in the hexagonal system: Acq is the 
equilibrium Lyapunov exponent for the closed system. The following quantities characterize the fractal repeller for 
viscosity with x — 60^/ri: Anoq is the nonequilibrium Lyapunov exponent of the repeller, /iks its KS entropy, 7 its 
escape rate, ci its partial information codimcnsion, and ch its partial Hausdorff codimcnsion. 



n 


Acq 


■^neq 


7 


ft-KS = Anoq — 7 


Cl = 7/ ^noq 


CH 


0.40 


1.5156 


1.5163 


0.0017 


1.5146 


0.0011 


0.0011 


0.50 


2.3519 


2.3539 


0.0023 


2.3516 


0.00098 


0.00092 


0.60 


3.7258 


3.7249 


0.0015 


3.7234 


0.00040 
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FIG. 10: Typical trajectory which moves through the whole system but which has a Helfand moment that remains close to 
zero. This trajectory belongs to the repeller of viscosity but not to the one of diffusion. 
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FIG. 11: Comparison between two methods of calculating the shear viscosity coefficient if = rj^y^xy in the hexagonal geometry: 
the Einstein-Helfand formula (continuous line) and the escape-transport formula llljl l with x ~ 60-yn (dots). 



B. Hausdorff dimension and viscosity 



In order to determine the viscosity by the chaos-transport formula (|25|l . we need to determine the partial Hausdorff 
codimension ch of the fractal repeller. The corresponding dimension du = 1 — ch is the Hausdorff dimension of the 
vertical asymptotes of the escape-time function depicted in Fig. [S] Its values range in the interval < dn < 1. 

The Hausdorff codimension can be obtained using the following numerical algorithm developed by the group of 
Maryland 0|. We consider an ensemble of pairs of trajectories starting from initial conditions (j)Q differing in a value 
e. The time taken by the trajectories to escape out of the system is given by the escape-time function in Fig. 
The pair is said to be uncertain if the trajectories and their Helfand moment present at least one of the following 
conditions: (i) both trajectories follow paths that differ by the successive passages through the cell boundaries, that is, 
if we associate to each trajectory a symbolic sequence (uji,lu2, ■■■) which gives the labels of the cell boundaries across 
which the successive passages occur, and both sequences are different; (ii) one of both trajectories has its Helfand 
moment which reaches the fixed absorbing boundaries H28|l when the Helfand moment of the other one still remains 
within these limits. If the pair does not present one of these conditions it is called certain. The fraction /(e) of 
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FIG. 12: Comparison between two methods of calculating the shear viscosity coefficient rj^y^xy in the square geometry: the 
Einstein-Helfand formula (continuous line) and the escape-transport formula 11311 with x ~ 45^71 for density n < 0.66, 
X = lOO^/n for 0.67 < n < 0.75, and x = 15Q^/n for 0.76 < n (dots). 




Q h I I I I I r 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 

density 



FIG. 13: Equilibrium Lyapunov exponent versus density in the hexagonal geometry. 



uncertain pairs in the initial ensemble is known to behave as the power 

f{e) ~ , (31) 

giving the HausdorfF codimension as its exponent. Derivations of this result can be found elsewhere 0,0,|23- This 
method has been already used in Refs. [iSllliL [l5l | . 

We have here applied the Maryland algorithm to obtain the HausdorfF codimension of the fractal repeller of viscosity. 
Table I compares the partial Hausdorff codimension with the partial information dimension in particular cases. We 
observe that both codimensions take very close values as expected. 

By varying the parameter we have obtained the shear viscosity thanks to the chaos-transport formula l|25|l. 
These values are plotted in Figs. 1151 and [TBI for the hexagonal and square geometries, respectively. We consider the 
shear viscosity coefficient rj* in the hexagonal geometry and the element ri,j.y xy* of the viscosity tensor in the square 
geometry. The values obtained with the chaos-transport formula (|25|l are compared with the values obtained by the 
escape-transport formula H13|l and those by the Einstein-Helfand formula Q obtained in the previous paper 17|. The 
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FIG. 14: Equilibrium Lyapunov exponent versus density in the square geometry. 



agreement between the three formulas is excellent, which confirms the theoretical results. 
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FIG. 15: Comparison between the three methods calculating the shear viscosity coefficient ry* in the hexagonal geometry: the 
Einstein-Helfand formula I^J (continuous line), the escape-transport formula H13|l (dots), and the chaos-transport formula H25|l 
(crosses) with x ~ 60^71. 



V. CONCLUSIONS 



In the present paper, we have applied the escape-rate formahsm to the computation of shear viscosity in the two-disk 
model by Bunimovich and Spohn jig . 

The escape-rate formalism implies the appearance of a fractal repeller associated with viscosity. We have numerically 
generated the fractal repeller associated with viscosity in this model and we have provided evidence for its fractal 
character. Using the chaos-transport formula of the escape-rate formalism, we have been able to evaluate the shear 
viscosity from the positive maximum Lyapunov exponent and the Hausdorff codimension of the fractal repeller of 
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FIG. 16: Comparison between the three methods calculating r]xy,xy in the square geometry: the Einstein-Helfand formula @ 
(continuous line), the escape-transport formula 1131 (dots), and the chaos-transport formula II25II (crosses) with x ~ 'i5\/n for 
density n < 0.66, x = lOOTn for 0.67 < n < 0.75, and x = 150^/n for 0.76 < n. 

viscosity. The values obtained by using the chaos-transport formula for shear viscosity have been compared with 
the values obtained by other methods based on the Einstein-Helfand formula, which is equivalent to the Green- 
Kubo formula as shown in our previous paper [l7j | . An excellent agreement has been observed between the different 
methods. This agreement brings an important support to the escape-rate formalism as a method to establish a 
connection between the transport properties - here of viscosity - and the underlying microscopic chaotic dynamics. 
The agreement therefore confirms the theoretical results of the escape-rate formalism 0, . 

This confirmation opens perspectives for our understanding of the hydrodynamic modes of shear viscosity. Indeed, 
in the case of diffusion, the fractal character of the repeller has been shown to imply that the nonequilibrium states 
are described by singular measures defined in the phase space of the system. These results suggest a similar scenario 
in the case of viscosity that the nonequilibrium states of a sheared fluid should also be described by singular measures 
at the phase-space level of description. 

Our results can be extended to investigate the connection between viscosity and underlying chaos in many-particle 
systems with a whole spectrum of positive Lyapunov exponents. In many-particle systems, the fractal repeller is 
characterized by a spectrum of partial fractal dimensions which enter the chaos-transport formula (|23|l . In this way, 
we can decompose a transport property such as viscosity onto the spectrum of Lyapunov exponents. 

The escape-rate formalism also opens perspectives to study viscosity and other transport properties in nanoscopic 
systems such as nanopipes or atomic and molecular clusters. Indeed, the escape-rate formalism provides a way to 
define the transport properties already in nanoscopic systems containing a very small number of particles. The escape- 
rate formalism is particularly appropriate for nanoscopic systems because the transport coefficients can be defined 
with the escape rate of the Helfand moment out of a finite range of variation. Since the Helfand moment for viscosity 
can be seen as the center of momenta of the particles its variations are limited to a finite range in a system containing a 
finite number of particles. In this way, the transport properties could be studied for nanoscopic systems of increasing 
size in a similar way as the equilibrium thermodynamic properties have been studied in nanoscopic systems as a 
function of their size. 

We hope to report on these issues in future publications. 
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